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ABSTRACT 



We provide simple polynomial fits to the X-ray photoelectric cross-sections 
(0.03 < E < lOkeV) for mixtures of gas and dust found in protoplanetary disks. 
Using the solar elemental abundances of Asplund et al. (2009) we treat the gas 
and dust components separately, facilitating the further exploration evolutionary 
processes such as grain settling and gain growth. We find that blanketing due 
to advanced grain-growth {ttmax > Ifim) can reduce the X-ray opacity of dust 
appreciably at Ex ~ IkeV, coincident with the peak of typical T Tauri X-ray 
spectra. However, the reduction of dust opacity by dust settling, which is known 
to occur in protoplanetary disks, is probably a more significant effect. The ab- 
sorption of 1-lOkeV X-rays is dominated by gas opacity once the dust abundance 
has been reduced to about 1% of its diffuse interstellar value. The gas disk estab- 
lishes a floor to the opacity at which point X-ray transport becomes insensitive 
to further dust evolution. Our choice of fitting function follows that of Morrison 
& McCammon (1983), providing a degree of backward-compatibility. 

1. Introduction 



The evaluation o 



(e.g. 



Brown &: Gould 



interstellar X-ray cross-sections stretches back over four decades 



19701) and is essential for interpreting X-ray observations subject 



to absorption by intervening material. Today, the evaluation and manipulation of highly 
detailed X-ray photoelectric cross-sections for arbitrary elemental abundances is made 
possible by powerful publicly available codeqj. Of some historical importance is the 



^These codes include XSPEC (Arnaud 



SHERPA flFreeman et al. 



2001 



1996 



http://heasarc.nasa.gov/xanadu/xspec/ ) , 



http://cxc.harvard.edu/sherpa/), and ISIS 
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paper by iMorrison &: McCaminonl (119831 . hereafter MM83) in which the authors present 
the photoelectric cross-section using a simple piece-wise quadratic polynomial fitted 
over discrete energy intervals (spanning a total range 0.03 < E < lOkeV). The main 
benefit of using a fitting formula is that it is an expeditious way to incorporate the 
cross-sections into a computer code. The basic limitation of this approach is that elemental 



abundances become locked in. This was relaxec 



somewhat i n the follow-up paper by 



Wilms et al 



teood hereafter WOO) made 



Balucinska- Church fc McCammonI ( 119921 ). while 

adjustments for the H2 molecule and included more accurate estimates of elemental 

abundances. Both these works treated the individual elemental abundances as adjustable 

parameters, increasing the versatility of the results whilst increasing the complexity of the 

implementation. 

Central to the piece-wise construction of MM83 is the fact that - despite their low 
cosmic abundances relative to H and He - the metals (C, O, Ne, Mg, etc.) contribute 
significantly to the X-ray opacity at energies above their respective K-shell photoelectric 
thresholds (0.28keV for C, 0.53keV for O, 0.87keV for Ne, 1.3keV for Mg, for example). 
These are the elements which are most readily incorporated into dust grains, although 
under typical interstellar conditions significant fractions of C and O (mostly in the form of 
CO), and Noble elements, persist in the gas phase. As such, the deposition of X-ray energy 
{E > 0.3keV) occurs through absorption by both gas and dust. The relative importance 
of gas and dust is highly energy dependent. At low energies, E < IkeV, only a handful 
of metals have sufficiently low K-shell thresholds that they can add significantly to the 
opacity fioor set by H+He. In this case most E ~ IkeV X-rays will be processed by the 
gas. While a greater variety of elements are susceptible to K-shell ionization at higher 
energies {E ~ lOkeV), the strong ~ E^^ dependence of the photoelectric cross-section 



(JHouck fc Denicola 



2000[ http://space.mit.edu/CXC/isis/) 
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means that the total cross-section drops rapidly with increasing energy. Here heavy metals 
such as Fe dominate the photoelectric cross-section, implying that lOkeV X-rays will likely 
be processed by grains. 

It is well established observationally that young pre-main sequence (T Tauri) stars are 
X-ray luminous with steady X-ray spectra typically peaking at energies of -E ~ 1 — 2keV 



( JFeigelson et al 



20051). In addition, there is considerable variability due to the eruption 
of stellar flares which generate a relatively hard X-ray component, briefly extending the 
X-ray spectrum to beyond lOkeV. In order to understand how X-rays interact with the disk 
it is necessary to consider how the evolution of protoplanetary disk material might affect 
the total gas-|-dust photoelectric cross-section. In this paper we extend the approach of 
MM83 by considering the effect on the X-ray photoelectric cross-section induced by two 
fundamental grain processes believed to occur in protoplanetary disks; dust settling and 
grain-growth. 

Interstellar dust grains become optically thick to X-rays {E ~ IkeV) once they 
grow to radii of a > O.Syum. At this point the interior atoms will be shielded from the 
impinging X-rays by the layers of atoms closer to the grain surface: the interior atoms do 
not see the full radiation eld and so do not contribute equa lly to the X-ray cross-section. 



This effect is referred to as self-blanketing (jFireman 



1974 WOO), rendering dust less 



efficient at absorbing X-rays. This reduction in absorbing efficiency is encapsulated in the 
self-blanketing parameter fb{E). As grains continue to grow to larger sizes (the precise 
size is highly energy dependent) essentially all X-rays incident upon them are absorbed 
(the regime of ge ometric optics). Ho wever, even the largest dust grains in the diffuse ISM 



0.25mm, 



Mathis et al. 



19771 hereafter MRN) are only marginally thick to IkeV 



X-rays, whereas the inferred grain-si ze distribution in protoplanetary disks can extend to 



signicantly larger sizes {amax > Imm, 



Throop et al 



2001 



Wilner et al 



2005; 



Lommen et al. 
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20091 ). making them strongly self-blanketing in the E = 1 — lOkeV range. Whether large 
grains actually contribute to the absorption of stellar X-rays depends largely upon whether 
they form part of the gas-dust mixture in the X-ray irradiated surface layers. 

In a steady protoplanetary disk the vertical component of the stellar gravitational eld, 
combined with the differential drift between the ballistic gr ains and partially centripetally 



supported gas, generates friction between the gas and dust (IWeidenschilling fc Cuzzi 



19931). 



The subsequent drag forces lead to vertical settling of dust toward the disk midplane, 
occurring on a timescale of order lO^years. This is short compared to typical disk lifetimes 
(> 10^ yrs ) , suggesting that disks exhib it a dust-poor atmosphere o verlying a thin, dust-ric h 



midplane (JDuUemond fc Dominik 



D'Alessio et al. 



(120061) 



20041 ). Following the notation of 
the extent of grain removal from the upper disk layers is represented by the parameter e, 
which we simply define as the local dust-gas mass ratio relative to that in the ISM (i.e. in 
the ISM e = 1). While dust settlin g occurs aggressive ly in typical T Tauri systems (the 



median value in Taurus is e ~ 0.01, 



Fur Ian et al. 



20061 ). the modeling of infrared spectral 



energy distributions suggests that there is a continual replenishment and recirculation of 



small grains in the upper layers of the disk (JDullemond fc Dominik 



20081 ). 



In this paper we return to the approach adopted by MM83, providing polynomial fits to 
the X-ray cr oss-sections usi n g the most recent estimates for the solar elemental abundances 



provided by 



Asplund et al. 



( 120091 ). By explicitly sphtting the total X-ray cross-section 
into gas and dust components it is possible to explore the effects of dust-specie processes 
(grain growth and settling). While the potential effect of grain-growth was mentioned in 
the aforementioned papers, it was shown to be of little quantitative importance for grains 
in the diffuse interstellar medium (typically affecting the cross-sections at the few percent 
level). In contrast, these effects are expected to become important in protostellar and 
protoplanetary systems where the nature of the gas-dust mixture departs signicantly from 
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that found in the interstellar medium ( iBeckwith et a 



mode 



1999 



ing of X-ray transp ort and energy deposition (JMaloney et al. 



2000f). Such results may facilitate the 



Nomura et al 



1996 : 



Igea fc Glassgold 



20071 ): resolving the various roles of X-rays in relation to other sources 



of en e rgy ( e .g. cosmic-rays an d ultraviolet irradiation 



2008 



2009 



Owen et al. 



Glassgold et al. 



2004 



Ercolano et al. 



20101): and 



sample disk material (iGiidel et al. 



interpreting the observed absorption of X-rays that 



2008 



Kastner et al. 



20051 ). We present separate 



polynomial fits for three components; 

1. A gas composed of hydrogen and helium. 

2. A gas composed of hydrogen, helium, the noble gases, and some fraction of the carbon 
and oxygen (required to make gas-phase CO). 

3. Dust grains containing all the elements heavier than helium, except the noble gases 
and a fraction of the carbon and oxygen. 



We regard Component 1 as a basal reference appropriate for chemically depleted gas. 
Component 2 is more typical of the gas in classical T Tauri disks. Component 3 is the 
dust, and it is to this that we apply our grain-specific effects. Components 2 and 3 together 
constitute disk material. The remainder of this paper is structured as follows. In Section 2 
we describe the fitting formula and the modifications required to include the grain-specie 
physics. The composition of the separate gas and dust components are given in terms of the 
latest solar abundances. In Section 3 the resulting X-ray cross-section fits for the various 
components are provided both in graphical and tabulated forms. The paper concludes with 
Section 4, providing a summary of the results. 
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2. X-ray cross-section fits 

In this section we explicitly separate the gas and dust components and provide 
polynomial fits to their cross-sections. To enable backward compatibility with MM83 
we have adopted the same fitting function and energy ranges. The X-ray photoelectric 
cross-section per H nucleus, <j{E), is described by the piecewise polynomial fitting function, 

a{E) = 10-2^ X (co + ciE + ca^^) E'^cm^. (1) 

where E is the X-ray energy and cq, ci, and C2 are the coefficients to be found. We exte nd 



Equa tion [T] by decomposing the medium into a mixture of gas and dust (e.g. 



Fireman 



19741 ). Initially we might split the total cross-section into two contributions so that 



ctot = o'gas + o'dust (units cm per H nucleus). From the point of view of computing X-ray 
opacities, the relatively small dust grains in the diffuse ISM can be treated as though their 
constituent atoms are in the gas phase (WOO). In this case the distinction between gas and 
dust can be suspended. As discussed in Section [1] the dust in protoplanetary disks may 
undergo both settling and growth. To accommodate both these phenomena we introduce 
two quantities, e and fb{E), encapsulating dust settling (the physical removal of dust) and 
grain growth (self-blanketing) respectively. These are grain-specific parameters, and as such 
only affect the contribution due to dust. Therefore we write the total cross-section as, 

o"tot = CTgas + efb{E)adust- (2) 

This is also true for combining fitting coefficients, 

cf = cr + eME)cr'\ (3) 



2.1. Elemental abundances and composition 



The elemental a bundances we have adopted are given in Table [T] (based on data from 



Asplund et al 



20091 ) ■ It is from these abundances that we make up our total gas-dust 
mixture. A comparison with previous estimates of the solar elemental abundances is 
given in Figure [H Although all the elements listed appear in the X-ray cross-section to 
some extent, carbon and oxygen (important contributors at -E ~ I keV) both have sola r 



abundance measurements that have varied considerably over time (lAsplund et al. 



20051). 



The distribution of elements between dust and gas phases is treated in a largely binary 
fashion; the metals are entirely incorporated into dust (/3 = 1), whereas H,He and the 
noble gases remain entirely in the gas phase (/3 = 0). The inclusion of Noble element s 



in the gas phase is rnotivat ed not only by observations and modeling (jPascucci et al.l 



2007 



Glassgold et al. 



20071 ). but also by their known low chemical 



reactivity and low 



200ll ). The exceptions to 



propensity for freezing-out onto grain mantles (iCharnley et al. 

the simple partitioning of elements ar e carbon an d oxygen, which are divided between 



the gas and solid phases according to 



Sofial fl2004h . Some C and O are 



necessary in the 



1991 : 



Thi et al. 



gas p hase to account for observations of circumstellar CO (iCalvet et al.l 

200lh . The partitioning of metals between gas and dust phases is an interesting avenue for 

furt her exploration, however we fix it in this paper to limit the number of free parameters 



fc.f. 



Balucinska- Church fc McCammon 



19921 ). The total cross-section of a mixture of 



elements (cm^ per H nucleus) is found by combinin g the abundanc es with the atomic X-ray 



cross-sections oz from the online NIST database g (iChantler 



I 



2000), 



o = T^z<^zA2 



(4) 



^http://www. nist.gov/physlab/data 
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In this paper we only consider neutral species, requiring that the ionization state of the 
gas is low. In order to produce simple fits to the cross-section we neglect solid state features 
near photoelectric thresholds. The detailed study of X-ray fine structure may p rovide a 



20091 ). As noted 



useful diagnostic of the composition of solid interstellar material (iLee et al.l 
in WOO, H2 has a cross-section that is approximately 40% larger (per H nucleus) than 
atomic H. We adopt the enhanced H2 cross-section and assume the gas is fully molecular 
(i.e. n{}i2)/n-^ = 0.5). Relative to earlier work focusing on the largely atomic ISM, the 
inclusion of the molecular enhancement increases the cross-section most signicantly at 
extreme ultraviolet wavelengths where hydrogen opacity dominates. 



3. Results 

The well-tested IDL tting routines POLY_FIT and SVD were used to evaluate and test 
the uniqueness of the ts. The fitting coefficients are given in Table [TJ The resulting opacity 
curves are shown in Figure [2l By comparison with a gas of pure H and He, Figure |2] clearly 
shows that C, O and Noble gases Ne and Ar, contribute signicantly to the gas opacity 
at energies E > 0.3keV. Away from photoelectric edges the fitting errors are typically 
at the 1% level, which is small compared to the compositional uncertainties of gas and 
dust. Noticeable in the dust component are significant (~ 10%) error spikes near to metal 
K-edge thresholds. These errors (consistent with those in WOO) arise from the numerical 
discretization of the elemental cross-sections, and from the fundamental limitations of a 
low-order fit. It is important to note that striving for a more accurate polynomial fit at 
K-shell edges is not entirely meaningful unless we include actual solid- state effects, which 



in general are too complex for a low order polynomial fit (JDraine 



2003h . 
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3.1. Grain growth, fb{E) 

The derivation of the self-blanketing factor fb{E) for spherical, homogeneous grains is 
described in Appendix A. Plots of fb{E) versus X-ray energy for a range of grain sizes are 
shown in Figure [31 As a function of energy it bears the imprint of cxdust since fb{E) depends 
upon the optical depth of the grain. At low energies {E <^ IkeV) the self-blanketing effect 
can be very strong (the grains are optically thick, fbiE) ^1), but it is not seen in the total 
(gas+dust) opacity because at these energies the total opacity is dominated by gas, not 
dust. At energies E > IkeV the metals in dust grains dominate the total (gas+dust) opacity 
but tend to be optically thin (for a < lOO/xm), and as a result the self-blanketing effect is 
weak. At intermediate energies {E ~ IkeV) the self-blanketing factor departs from unity 
once grains have grown to a ~ 0.5/im. Grains as large as a ~ 1mm are extremely optically 
thick, hiding 99.9% of their constituent atoms from the ambient X-rays. The grain-size 
distribution in protoplanetary disks is largely unknown, although the process of coagulation 
will tend to drive mass into larger grains. The ensemble averaged self-blanketing factor, 
< fb{E) >, for an population of grains with a power-law size distribution, dn/da oc a~^, is 
given by the average, 

< ME) >= JaLJ_,^^ • (5) 

where p = 3.5, Omm ~ O.Olfim and amax ~ 0.25/im for the classic MRN distribution. 
The mass per decade of grain size in the MRN distribution is weakly dominated by the 
largest grains. This becomes increasingly true if we flatten the distribution through grain 
coagulation (taking mass from the small grains and placing it in the big grains). It is 
therefore of interest to briefly study the effects of changing the maximum grain size amax 
and power-law index p. In Figure H] we preserve the slope of the MRN distribution but 
allude to grain-growth by gradually increasing amax from the interstellar value ~ 0.25/xm 
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to 1mm. The increase in amax draws an increasing large fraction of the dust mass into the 
realm of self-blanketing. A similar effect can be achieved by reducing the power-law index 
p, as shown in Figure O In this case the evaluation of < fb{E) > is increasingly dominated 
by grains with sizes a ~ CLmax (we hold amax = 10/xm so that the curves for p < 2 can be 
compared to the curve for the single grain size population a = lO/im in Figure [3]). Fitting 
coefficients for the dust component computed using MRN distributions {amax = 1, 10, 100 
and 1000/im) are given in Table [2j 

The effect of self-blanketing on the total (gas+dust) photoelectric cross-section is 
shown in the bottom panel of Figure [Hi The onset of grain-growth initially affects the 
cross-section near Ex ~ IkeV - coincident with the peak X-ray emission from T Tauri 
systems. 

3.2. Dust settling, e 

The effect of dust settling is shown in the upper panel of Figure [6l and is straightforward 
to interpret. To isolate this effect we show its affect applied to the standard interstellar 
MRN grain size distribution, for which self-blanketing is mild due to the preponderance 
of small grains (a < 0.5/im). Once e has dropped to below 0.1 (consistent with T Tauri 
systems in Taurus) the contribution made by dust to the total opacity at E = IkeV is 
small compared to that of the gas. The gas sets a basal floor to the cross-section. Further 
reduction of the cross-section is only possible by depleting gas phase metals and noble 
elements. 
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3.3. X-ray scattering 



The primary purpose of this paper is to quantify the photoelectric absorption 
cross-sections that determine how X-ray energy is deposited into gas and dust. The 
more general consideration of how X-rays are transported through space requires the 
inclusion of scattering. Restricting the discussion to scattering by spheres of radius 
a, the main parameters that determine scattering efficiencies are the size parameter 
X = 2TTa/ X where A is the X-r ay wavelength, and the real part of the complex refractive 
index, m. Ivan de Hulstl ( 119571 ) provides an authoritative discussion of scattering processes 
throughout this parameter space (see Fig. 20 Ch. 10). While Mie theory provides a 
rigorous way to compute the extinction of individual particles, a number of approximations 
are applicable at the limits of this parameter space. The region of interest is reduced 
somewhat since m — 1 ^ 1 for r nost materials (including interstellar dust candidates) at 



X-ray wavelengths {E > O.lkeV, 



Draine 



20031 ). A refractive index close to unity implies 



that the phase of an X-ray inside the grain only slowly deviates from that of an exterior 
X-ray. When the total differential phase change 2x{m — 1) is smaller than unity (true for 
most interstellar grai ns) one may use th e Rayleigh-Gans approximation to compute the 



scattering efficiency (JMathis fc Lee 



199lh . The Rayleigh-Gans approxir nation asserts th at 



19651 1. 



the grain can be viewed as an ensemble of Rayleigh scattering centers (lOverbecld 

The sum total of these fields gives a scattered field that is sharply peaked in the forward 

directio n: the property responsible for the narrow halos observe d around X-ray point 



sources ( 



Roll 



1983 



Predehl fc Schmitt 



1995 



Smith fc Dwek 



19981 ). When computed for the 



^"""""" '' ""'"^ ^ '"'" P0P«lation, 90% of IkeV X-.ay. „e scattered by le. 
than one degree away from the forward directioio. 



■^In the limit of geometrical optics (large grain sizes, 2x(m — 1) ^ 1) the absorption and 
scattering (diffraction) cross-sections are each equal to the geometrical cross-section vra^ (c.f. 
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In Figure [7] we show alongside the photoelec tric absorption cross-section the scattering 



Draind (120031 ) . When the E^ factor has 



cross-sections of both gas and dust computed by 

been taken into account it can be seen that the non-relativistic (Thomson) scattering 

by gas is essentially wavelength independent. Thomson scattering - the interaction of 



a photon with a free (or weakly bound) electron - is due alr nost entirely to t 



le H and 



He since these elements carry the majority of the electrons (Jgea fc Glassgold 



19991 ). It 



only contributes significantly to the total extinction ai E > 6keV. Even in dust depleted 
scenarios, scattering by gas is only ever important for the highest energy photons observed 
in T Tauri spectra. In contrast, dust scattering contributes significantly at lower energies, 
E ~ IkeV. However, while small-angle scattering by grains is important for interpreting 
observations of X-ray halos, the very narrow phase function associated with dust scattering 
suggests that multiple scattering events are required to produce a significant diffuse X-ray 
field inside a protoplanetary disk. Thus, while the dust scattering cross-section may be 
large in magnitude it does not redirect radiation effectively. Since the total photo-electric 
absorption efficiency at -E ~ 1 — 2keV is never much smaller than the scattering efficiency, 
the majority of X-ray photons emitted by T Tauri stars will have been absorbed before a 
significant diffuse component is generated. 

4. Summary 

We have used the most recent estimates of the solar abundance to construct simple 
parametric fitting functions for the X-ray photoelectric cross-sections of material in T Tauri 
disks. By separating the medium into gas and dust components we have shown how the 
X-ray opacity of T Tauri disk material will change as dust grows in size or settles towards 



extinction paradox). 
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the disk midplane. Here we provide a summary of our results: 

1. Separate polynomial ts are presented for the photoelectric absorption cross-section 
of gas and dust (Tables [1] and [21 Fig. [2]). These fits are backward-compatible with MM83. 

2. Grains as large as a ~ 1/im are required before the effect of self-blanketing 
is seen. This effect is most readily seen at -E ~ 1 — 2keV (Fig. [3]) - coinciding with 
the peak X-ray emission of T Tauri systems. Due to the large fraction of small grains 
(a << 1/im), considerable grain-growth is required {amax » lAini) before the MRN 
grain-size distribution exhibits appreciable self-blanketing at -E ~ IkeV (Fig. HI). 

3. Provided sufficiently large grains are present, the self-blanketing factor averaged 
over a power-law grain-size distribution is sensitive to the power-law index (Fig. \5^. For 
size distributions flatter than MRN (i.e. p < 3.5) we observe a stronger ensemble-averaged 
self-blanketing effect, eventually asymptoting to that of a single-sized grain population with 

4. Realistic structural and compositional inhomogoneities of grains will change the 
self-blanketing results. Non-spherical grains will self-blanket less effectively. The layering 
of grain mantles may either increase or decrease the self-blanketing factor, depending upon 
the composition and order of these layers. 

5. The degree of dust settling evident in protoplanetary disks, e ~ 0.01, is more than 
sufficient to reduce the X-ray opacity of dust to levels that are small compared to those of 
the gas (c.f. Fig. |6]). The combined grain-growth and settling observed in protoplanetary 
disks suggests that X-ray energy in T Tauri systems is preferentially injected into the gas 
component. 

6. While both grain-growth and settling reduce the X-ray cross- sect ion, the remaining 
gas opacity sets a lower limit that ensures atot {E ~ IkeV) does not drop below about 50% 
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of the diffuse ISM value (Fig. E>- 

7. Further reduction of the X-ray cross-section in dust-evolved systems is contingent 
upon changes in the metal composition of the gas. The removal of CO and Ne as the gas 
disk dissipates could potentially reduce atot by a further factor of three. 

X-ray counterparts have been identified for many protoplanetary disk systems in Orion, 
several of w hich display X-ray spectra consistent wit h absorption by surrounding material 



(e.g. COUP 



Getman et al. 



2005 



Kastner et al 



20051). By comparing the absorption column 



inferred from X-rays to that determined from visual extinction (by dust), it is in principle 
possible to place constraints on the composition and dust-to-gas ratio of intervening disk 
material. Edge-on sources COUP 419 and COUP 241 are such examples, suggesting that 
the absorbing material is considerably deficient in dust (an order of magnitude or more 
relative to to interstellar value). 



Giidel et al. 



(120081 ) have detected bipolar X-ray jets emanating from the strongly 
accreting T Tauri system DC Tau. Comparison of the two jets indicates that one is probably 
viewed through an excess column of material {Nh ~ 2.7 x lO^^cm"^) due to the surrounding 
disk. Unfortunately, while this scenario potentially affords a differential analysis of the 
intervening material, the visual extinction along this line of sight is uncertain, rendering the 
gas-to-dust ratio determination inconclusive. 

Far-ultraviolet radiation (FUV) is another important source of energy in protoplanetary 
disks. The transport of continuum-FUV is largely controlled by dust, and as such is highly 
sensitive to dust evolution. In particular, as dust settles towards the disk midplane the 
upper parts of the disk become increasingly transparent to FUV photons, reducing the FUV 
opacity to X-rays levels. In contrast, the propagation of X-rays is impeded almost entirely 
by gas after only a modest degree of dust evolution (e ~ 0.1). The deposition of X-rays with 
energies typical of the steady coronal emission from T Tauri stars (E' ~ 1 — 2keV) will be 



de posited in the uppe r 
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layer s of the disk - the so-called warm molecular layer. In the models 



of lAikawa fc Nomura! (l2006l ) the warm layer of gaseous CO at R=200AU has a typical 
vertical column density of A^^(CO)~ lO^^cm"^. We expect this layer to be very optically 
thick, Tikev ~ 100, when viewed at the small grazing angles (~ 5 degrees) appropriate for 
impinging stellar X-rays. The impenetrability of this layer is mitigated somewhat at higher 
energies {E ~ lOkeV) where the photoelectric cross-sections becomes small compared to 
the scattering cross-section . In this case, Thomso n scattering by gas can be an effective 



means of disk penetration f jlgea fc Glassgold 



19991 ). 



TB and EB gratefully acknowledge funding by NASA under grant NN08 AH23G from 
the ATFP and SSO programs. 

A. Derivation of dust self-blanketing factor fbiE) 

A single interstellar dust grain of sufficient size may become optically thick to X-rays. 
In such a case the interior atoms see a reduced flux of X-rays relative to the atoms on 
the surface of the grain. The total absorptive efficacy of the dust grain is therefore less 
than that of an equal mass of compositionally identical gas. This effect is referred to as 
self-blanketing, and to calculate its magnitude we must consider the dust grain as a radiative 
transfer problem by itself (albeit a highly simplified one). In this Appendix we consider 
the self-blanketing of spherical dust grains - not because the sphere is a physically sensible 
model for interstellar grains but because the sphere is the most self-blanketing geometry 
(under isotropic illumination). The self-blanketing effect of real (i.e. non-spherical) grains 
will be less than that presented here. We adopt the basic notation of WOO, although we 
differ slightly in our approach by explicitly specifying a grain shape (sphere). The average 
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cross-section of the constituent atoms in the sohd material is given by 

-(p\ _ ^zAzPz/^zjE) _ adust /..x 

The atomic abundances Az and depletion factors (3z are provided in Table 1. In effect 
Eqn Al simply removes the per hydrogen scaling, replacing it with per grain atom. The 
impingent intensity Iq is partially absorbed by the grain, emerging with an intensity I{x) 
depending upon the offset, x, between the ray and the grain center. If the grain has a 
number density n then the ray must traverse an optical depth 



TiX 



2anVa2 - x"^. (A2) 



Averaged over the geometric cross-section of the grain the average emergent intensity is 

< J > 2 /"" / \ 2 

— - — = ^ / xexp ( —2anVa^ — x"^ ) da; = ^ [1 — (tq + 1) exp(— tq)] (A3) 

h a Jo ^ ' ^0 

The final solution is expressed in terms of the diametrical optical depth, tq = 2ana. Since 

spheres look the same from all directions Eqn IA3I is already angle-averaged. The fraction of 

radiation removed from the beam is then simply 1— < / > //q. If we ignore self-blanketing 

and assume all the constituent atoms in the grain see the same incident flux of photons, 

the fraction of photons removed is f = 2/3ro. Comparing this with the self-blanketed result 

(Eqn. IA3p yields the self-blanketing factor for a spherical grain, 

3l-^[l-(ro + l)exp(-ro)] 
hm = ^ -^ (A4) 

2 To 

Note that a grain with a small self-blanketing effect has a large fb{E). In WOO the grain 
shape is not explicitly involved in the calculation. Rather, they characterize the grain solely 
by a notional optical depth, f and their self-blanketing factor is given as, 

r To 
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In the last step we have expressed their result using the appropriate mean optical depth 
for a sphere, f = 2/3ro. The difference between the WOO expression and our sphere-specie 
result (Eqn. IA4p is less than 5% in the range tq = 1 — 10. Both functions have the same 
asymptotic behaviour. For most purposes the computationally simpler WOO expression 
given by Eqn. IA5I will suffice. However, for grain-shapes less regular than the sphere the 
self-blanketing factor should be evaluated in a shape-specie manner. Grain structures with 
low volume illing factors will be less self-blanketing. 
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Fig. 1. — Comparison of the solar elemental abundances of 



Asplund et al. 



fl2009[ ) with those 



used previously in the literature. The quantity f{X) is the fractional abundance of e 



Asplund et al. 



ement 



(120051 ). 



X relative to hydrogen nuclei. Diamonds - abundances adopted by 

Plusses - abundances adopted by MM83. Of particular note are the variation s in the carbon 

and ox ygen abundances, which over time have varied by > 50% relative to the 



Asplund et al. 



(120091) values. 
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7ig. 2 . — Bottom - X-ray cross-section of our basic gas-dust mixture using 



Asplund et al. 



( 20091 ) solar abundances. The parametric fitting coefficients for these curves are given in 
Table 2. The ducial dust contribution is calculated without modication by grain processes 
(i.e. e = fb{E) = 1). To include the effects of settling and self-blanketing one must combine 
the gas and dust components following Eqns [2] and [31 The reader should note the E^ scaling 
- the total cross-section in fact drops with increasing energy. Top - Residuals of our fitting 



functions for gas and dust components. 
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Fig. 3. — Self-blanketing factor for spherical dust grains of radius a. At Ex = IkeV a dust 
grain with a ~ 0.5;um has an diametrical optical depth of r ~ 1, marking the onset of 
appreciable self-blanketing. 
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Fig. 4. — Self-blanketing factor averaged over a power-law size distribution of grains. The 
different lines correspond to changing the maximum grain size in the population, amax , 
while keeping the minimum size fixed, amin = 0.01/im. The power-law index of p = 3.5 and 
ttmax = 0.25/im is consistent with the classic MRN result for diffuse interstellar dust. 
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Fig. 5. — Self-blanketing factor averaged over a power-law MRN size distribution of grains. 
The different lines correspond to changing the power-law index of the population, p, while 
keeping the minimum and maximum grain sizes fixed at amin = O.Olyum and a^ax = lO/xm 
respectively. The power-law index of p = 3.5 is the classic MRN result for diffuse interstellar 
dust. We have extended amax to a value larger than that specied by MRN, since some grain 
growth is required before self-blanketing manifests itself. 
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Fig. 6. — Effects of grain processes on the total (gas+dust) X-ray photoelectric cross-section. 
Top - the contribution to the X-ray opacity made by grains varies in proportion to changes 
in the dust:gas mass ratio, e, occurring as a result of dust settling. Bottom - grain growth 
affecting the dust opacity in a more complicated, energy-dependent manner. In both plots 
the gas opacity sets a basal limit. 
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Fig. 7. — Comparison of the (gas+dust) photoelectric absorption {solid lines) and scattering 
due to gas {dashed line) and dust {dotted line), as a function of the setthng parameter e. 
At E r^ IkeV scattering by gas contributes neghgibly to the total cross-section, whereas 
dust scattering is of the same order of magnitude as dust absorption. At Ex ~ lOkeV the 
scattering opacity is large with both gas and dust contributing significantly. Although the 
gas may play a role in scattering X-rays at these energies, the total photoelectric absorption 
is dominated by the metals in dust. Note the E^ scaling - the (Thomson) gas scattering 
opacity is in fact almost independent of energy 
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Table 1. Elemental abundances and fraction in dust 

Element Abundance Dust Fraction^ 

12 + logAz Az = nz/nH Pz 

H 12.0 1 

He 10.93 8.51 x 10"^ 

Li 1.05 1.12 X lO'ii 1 

Be 1.38 2.40 x 10"" 1 

B 2.70 5.01 X 10"^° 1 

C 8.43 2.69 X 10"^ 0.28^ 

N 7.83 6.76 x 10"^ 1 

O 8.69 4.89 X 10"^ 0.50^ 

F 4.56 3.63 x 10"^ 1 

Ne 7.93 8.51 x 10"^ 

Na 6.24 1.74 x lO"*^ 1 

Mg 7.60 3.98 X 10"^ 1 

Al 6.45 2.81 X 10^"^ 1 

Si 7.51 3.24 X 10"^ 1 

P 5.41 2.57 X 10"^ 1 

S 7.12 1.32 X 10-5 1 

CI 5.50 3.16 X 10-^ 1 

Ar 6.40 2.51 x 10"^ 

K 5.03 1.07 X 10-^ 1 

Ca 6.34 2.19 x lO"*^ 1 
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Table 1 — Continued 



Element Abundance Dust Fraction^ 

12 + log A^ Az = nz/nH Pz 



Sc 


3.15 


1.41 X 


10-9 




Ti 


4.95 


8.91 X 


10-8 




V 


3.93 


8.51 X 


10-9 




Cr 


5.64 


4.37 X 


10-^ 




Mn 


5.43 


2.69 X 


10-7 




Fe 


7.50 


3.16 X 


10-5 




Co 


4.99 


9.77 X 


10-8 




Ni 


6.22 


1.66 X 


10-6 





^The fraction of element Z in grains. 



Sofial mm 
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Table 1. Fitting coefficients for different components 



Energy Range H2+He 

(keV) Co ci 



C2 



gas 



Co 



Cl 



C2 



dust 



Co 



Cl 



C2 



0.030-0.055 
0.055-0.100 
0.100-0.165 
0.165-0.284 
0.284-0.400 
0.400-0.532 
0.532-0.708 
0.708-0.867 
0.867-1.303 
1.303-1.840 
1.840-2.471 
2.471-3.210 
3.210-4.038 
4.038-7.111 
7.111-8.331 
8.331-10.00 



14.3 722 -4190 

22.3 432 -1530 
31.5 246 -578 

42 120 -198 

51.5 50 -68.7 

58.1 16.8 -26.4 

63.2 -2.37 -8.44 
66.8 -12.5 -1.18 

69.4 -17.5 1.2 

71.6 -22.3 3.56 
66.6 -16.8 2.06 
61.4 -12.6 1.21 
56.6 -9.6 0.75 

48.4 -5.8 0.308 

40.5 -3.43 0.128 
37.8 -2.8 0.091 



14.2 727 -4130 

22 445 -1550 

31 263 -614 

43.7 112 -165 

49 86 -103 

58.6 36.9 -39.9 

48 130 -82.2 

77.4 46.3 -22 

80.1 69.8 -28.3 

117 7.43 -1.87 

107 16 -3.75 

106 13.6 -2.63 

138 -1.99 -0.179 

142 -4.7 0.239 

138 -3.36 0.133 

88.9 8.15 -0.547 



0.0344 -1.62 88.2 

-0.147 4.19 48.1 

-0.677 14.9 9.6 

-1.12 23.6 -16.2 

0.188 24.6 -1.09 

-3.57 55.5 -37.9 

-8.24 89.6 -48.1 

57.1 -49.9 52.1 

9.11 72.7 -20.8 

-8.71 106 -25.7 

34.9 72.4 -11.4 

23.6 85.1 -11.3 

116 28.2 -2.55 

191 -2.92 1.09 

812 -74.7 6.49 

-33 137 -6.39 
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Note. — The coefficients belong to the fitting function given by Eqn [T] and can 
be combined according to Eqn OThe gas component consists of H, He, Ar, Ne and a 
fraction (1-/3^) of O and C (see Table 1). The dust component consists of the remaining 
elements. The H+Hc case is meant to serve as a basal reference. 



Table 2. Fitting coefficients for dust component subject to grain growth 



Energy Range 


amax = Iptm 


Ctmaa 


; = lO^rn 




(^raaa 


; = lOO^m 


(^max 


= lOOO^ni 




(keV) 


Co 


Cl 


C2 


C() 


Cl 


C2 


CO 


Cl 


C2 


CO 


Cl 


C2 


0.030-0.055 


0.0203 


-1.5 


47 


0.00785 


-0.576 


17 


0.00260 


-0.19 


5.52 


0.000829 


-0.0607 


1.75 


0.055-0.100 


0.0192 


-1.64 


50.9 


0.0166 


-0.918 


20.6 


0.00616 


-0.325 


6.9 


0.00203 


-0.106 


2.21 


0.100-0.165 


-0.173 


1.69 


40.7 


0.00492 


-0.794 


21.7 


0.00703 


-0.371 


7.65 


0.00275 


-0.129 


2.48 


0.165-0.284 


-1.36 


15.8 


2.1 


-0.236 


2.05 


14.4 


-0.0475 


0.28 


6.03 


-0.0124 


0.0528 


2.04 


0.284-0.400 


-2.96 


27.5 


-0.943 


-0.61 


4.29 


16.7 


-0.103 


0.539 


7.52 


-0.0238 


0.091 


2.58 


0.400-0.532 


-9.08 


61.6 


-37.5 


-3.24 


17.7 


2.89 


-0.755 


3.89 


4.26 


-0.212 


1.06 


1.66 


0.532-0.708 


-19.2 


100 


-48.7 


-7.71 


31.9 


0.587 


-1.8 


7.03 


4.67 


-0.506 


1.92 


1.9 


0.708-0.867 


17 


19.2 


17.5 


2.13 


8.01 


23.4 


1.92 


-2.27 


13 


0.725 


-1.17 


4.65 


0.867-1.303 


-14 


98.2 


-28.7 


-25.5 


71.1 


-11.5 


-5.4 


14.6 


3.75 


-1.4 


3.74 


1.93 


1.303-1.840 


-27.1 


120 


-28.6 


-63 


123 


-24.5 


-17.8 


31.3 


-0.219 


-4.69 


8.05 


1.04 


1.840-2.471 


15.8 


84.1 


-13.4 


-59.6 


119 


-17.9 


-28 


40.5 


-0.463 


-6.39 


9.35 


1.39 


2.471-3.210 


10.1 


91.3 


-12.1 


-63.7 


121 


-15.4 


-63.4 


66.4 


-4.23 


-14.7 


15.2 


0.66 


3.210-4.038 


106 


31.9 


-2.95 


40.9 


56.7 


-5.52 


-71.8 


72.3 


-5.27 


-15.4 


15.7 


0.556 


4.038-7.111 


187 


-1.84 


1.01 


156 


6.49 


0.416 


-1.51 


43 


-1.96 


-17 


16.9 


0.479 


7.111-8.331 


800 


-72.4 


6.37 


704 


-53.9 


5.39 


183 


37.1 


0.832 


24 


13.7 


2.21 


8.331-10.00 


-39.5 


138 


-6.42 


-90.7 


145 


-6.69 


-397 


181 


-7.76 


-276 


86 


-2.05 



00 



Note. — The coefficients apply only to the dust component, amax refers to the maximum grain-size in a MRN grain-size distribution. The 
total dust mass is held constant. 



